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Abstract 

Multi-scale  interactions  between  initially  planar  deformation  waves  in  heterogeneous  energetic  solids 
and  macro-scale  boundaries  were  computationally  examined  to  characterize  wave  interaction  struc¬ 
tures  and  dissipative  heating  responsible  for  combustion  initiation.  The  macro-scale  response  was 
described  by  a  continuum  theory  that  accounts  for  elastic  and  inelastic  volumetric  deformation  in  a 
thermodynamically  consistent  manner;  the  model  equations  were  numerically  integrated  in  terms  of 
generalized  curvilinear  coordinates  using  a  high-resolution  shock-capturing  technique.  The  meso- 
scale  scale  response  was  described  by  conservation  principles  and  an  elastic- viscoplastic  constitutive 
theory  to  predict  contact  induced  nonlinear  deformation  and  thermomechanical  fields  within  parti¬ 
cles;  the  model  equations  were  numerically  integrated  using  a  combined  finite  and  discrete  element 
technique.  Predicted  meso-scale  fields  were  locally  averaged  and  compared  to  independent  pre¬ 
dictions  given  by  the  macro-scale  model.  Wave-boundary  interactions  resulted  in  dispersed  wave 
structures  that  are  analogous  to  Mach  stems  in  gas  dynamics,  with  maximum  dissipative  heating 
rates  occurring  within  the  stem  region  along  bounding  surfaces  away  from  the  stagnation  region. 
These  maximum  heating  rates  substantially  increased  the  fraction  of  mass  locally  heated  to  ele¬ 
vated  temperature.  Bulk  shear  was  shown  to  have  a  small  effect  on  dissipative  heating  for  the 
impact  conditions  investigated  in  this  study.  The  dependence  of  wave  structure,  dissipative  heating 
rate,  and  hot-spot  mass  fraction  on  initial  wave  strength  was  characterized.  Macro-  and  meso-scale 
predictions  agreed  well,  provided  that  the  averaging  volume  size  was  suitably  chosen. 

1  Significance 

Interaction  between  deformation  waves  in  heterogeneous  energetic  solids  and  macro-scale  bound¬ 
aries  may  trigger  prompt  vigorous  combustion  due  to  locally  high  dissipative  heating  rates.  For 
example,  Wilson,  et  al  [11],  experimentally  showed  that  interaction  between  an  initially  planar  shock 
in  the  explosive  PBXW-128  (60%  HMX  and  40%  polymeric  binder)  and  a  hemispherical  glass  anvil 
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Figure  1:  Illustration  of  the  experimental  test  geometry  used  to  characterize  the  combined  shock- 
shear  behavior  of  explosives. 


initiated  vigorous  combustion  near  the  explosive-anvil  boundary  away  from  the  high-pressure  stag¬ 
nation  region;  an  illustration  of  the  test  configuration  is  shown  in  Fig.  1.  This  observation  suggests 
that  macro-scale  shear  may  be  important  for  combustion  initiation,  though  its  precise  role  remains 
unclear.  Analysis  of  deformation  wave-boundary  interactions  is  warranted  because  they  commonly 
occur  in  practice,  but  such  analysis  is  complicated  by  the  existence  of  multi-scale  processes  that 
arise  due  to  material  heterogeneity;  the  initial  distribution  of  particle  sizes  may  range  from  0.001- 
200  /im,  whereas  the  size  of  engineered  materials  is  orders  of  magnitude  larger.  Also,  length  and 
time  scales  of  physical  phenomena  occurring  at  the  particle  and  macro-scales  may  vastly  differ; 
dissipative  heating  can  occur  rapidly  over  length  scales  that  are  much  smaller  than  the  average 
particle  size,  whereas  detonation  of  explosives  resulting  from  low  pressure  impact  can  occur  slowly 
over  comparatively  large  distances  4  cm).  Though  it  is  well-accepted  that  material  heterogeneity 
gives  rise  to  deformation  wave  structures  exhibiting  intense  particle-scale  fluctuations,  the  relative 
importance  of  dissipative  mechanisms  and  their  influence  on  macro-scale  behavior  remain  largely 
unresolved  issues,  particularly  for  problems  involving  wave-boundary  interactions.  Consequently, 
multi-scale  computational  examination  of  such  interactions  is  warranted,  particularly  in  the  absence 
of  particle-scale  data. 

In  the  remaining  sections  of  this  report,  key  objectives  of  this  work,  and  research  methods  used, 
are  first  summarized,  followed  by  a  discussion  of  key  results.  A  list  of  research  personnel  supported 
by  this  work,  and  of  publications  and  presentations  which  resulted  from  it,  are  then  given. 

2  Objectives  and  Research  Method 

Motivated  by  the  experiments  of  Ref.  [11],  an  overall  objective  of  this  study  was  to  computationally 
examine  multi-scale  dissipative  heating  induced  by  the  interaction  of  initially  planar  deformation 
waves  in  granular  explosive  with  macro-scale  boundaries  and  to  characterize  its  dependence  on 
initial  wave  strength.  To  this  end,  computations  were  independently  performed  using  both  macro¬ 
scale  and  meso-scale  model  descriptions  to  obtain  a  better  understanding  of  this  process.  The 
macro-scale  description  can  be  applied  to  engineering  scale  systems,  but  does  not  explicitly  resolve 
particle-scale  temperature  fluctuations  (referred  to  as  hot-spots)  that  are  important  for  combustion 
initiation.  The  meso-scale  description  resolves  particle-scale  fluctuations,  but  cannot  be  faithfully 
applied  to  engineering  scale  systems  due  to  computational  limitations.  This  work  is  a  preliminary 
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step  towards  a  longer  term  goal  of  using  meso-scale  computations  to  facilitate  development  of  im¬ 
proved  macro-scale  constitutive  theories  for  impact  induced  phenomena  occurring  in  heterogeneous 
energetic  solids.  It  is  noted  that  the  meso-scale  modeling  effort,  which  represents  a  significant 
component  of  the  work  performed,  was  not  part  of  the  original  scope  of  the  proposed  project.  This 
supplemental  effort  was  undertaken  to  provide  information  about  deformation  particle-boundary 
interactions,  including  particle-scale  dissipative  heating,  in  the  absence  of  experiments.  Combined 
shock-shear  experiments  for  granular  explosive  planned  by  AFRL-MNME  personnel  that  were  to 
guide  development  of  a  macro-scale  constitutive  theory  for  this  project  were  suspended  due  to  pro¬ 
grammatic  issues.  The  objectives  and  research  method  for  both  the  macro-scale  and  meso-scale 
modeling  and  computational  efforts  are  summarized  below. 

2.1  Macro-Scale  Modeling  and  Computations 

Objectives.  Specific  objectives  of  the  macro-scale  modeling  and  computational  effort  included  the 
following: 

1.  Analysis  of  Deformation  Wave  Interactions  with  Planar  Boundaries.  (See  Ref.  [6]  for  a  com¬ 
prehensive  discussion.)  The  objective  of  this  study  was  to  computationally  examine  inter¬ 
action  between  an  incident  planar  deformation  wave  in  granular  HMX  explosive  (octahydro- 
l,3,5,7-tetranitro-l,3,5,7-tetrazocine)  and  a  rigid  planar  boundary  that  deflects  the  flow  through 
an  angle  6,  as  illustrated  in  Fig.  2(a).  Granular  HMX  was  chosen  because  its  uniaxial  com¬ 
paction  behavior  is  reasonably  well-characterized  [3,  4].  The  planar  boundary  represents  a 
theoretically  simpler  configuration  for  assessing  deformation  wave-boundary  interactions  than 
those  associated  with  more  complex  boundary  geometries.  Particular  emphasis  was  placed  on 
characterizing  how  incident  deformation  wave  strength,  expressed  in  terms  of  a  piston  impact 
speed  Up,  and  flow  turning  angle  6  affect  interaction  structure  and  dissipative  heating  of  the 
explosive,  where  0°  <  6  <  90°.  The  limiting  values  0  =  0°  and  90°  result  in  no  boundary 
interaction  and  pure  reflection,  respectively.  This  computational  analysis  also  enabled  the  ef¬ 
fects  of  rapid,  repeated  loading  to  be  investigated  as  material  pre-compacted  by  the  incident 
planar  deformation  wave  was  subsequently  deformed  by  reflected  waves.  As  such,  a  secondary 
objective  of  this  effort  was  to  examine  how  rapid,  repeated  loading  of  explosive  affects  its  sen¬ 
sitivity;  this  represents  a  largely  unanswered  question  that  has  important  ramifications  in 
practice. 

2.  Analysis  of  Deformation  Wave  Interactions  with  Non-Planar  Boundaries.  (See  Refs.  [5,  6] 
for  a  comprehensive  discussion.)  The  objective  of  this  study  was  to  computationally  examine 
interaction  between  an  incident  planar  deformation  wave  in  granular  HMX  and  a  rigid  semi¬ 
circular  boundary.  This  configuration  was  chosen  to  mimic  the  experimental  configuration 
of  Ref.  [11],  depicted  in  Fig.  1;  the  computational  domain  used  in  this  study  is  shown  in 
Fig.  2(b).  For  this  configuration,  the  flow  deflection  angle  decreases  with  position  along  the 
boundary  surface.  As  with  the  planar  boundary  study,  emphasis  was  placed  on  characterizing 
how  deformation  wave-boundary  interaction  structure  and  explosive  dissipative  heating  varies 
with  incident  wave  strength. 

Research  Method.  (A  comprehensive  discussion  of  the  model,  numerical  technique,  and  its  verifi¬ 
cation  and  validation  for  granular  HMX,  is  given  in  Refs.  [3,  4,  6].)  The  macro-scale  model  used  in 
this  study  is  an  extension  of  a  continuum  field  theory  for  the  Deflagration-to-Detonation  Transition 
(DDT)  of  granular  explosive,  first  formulated  by  Baer  and  Nunziato  [1],  and  later  improved  by 
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Figure  2:  Computational  domains  used  to  investigate  the  interaction  of  initially  planar  macro-scale 
deformation  waves  in  granular  explosive  with  (a)  rigid  planar  boundaries  and  (b)  rigid  circular 
boundaries. 


Bdzil,  et  al  [2],  to  account  for  both  elastic  and  inelastic  volumetric  deformation  in  a  thermody¬ 
namically  consistent  manner  [3,  4].  For  inert  impact,  e  =  pg/ ps  ^  1,  where  pg  and  ps  are  the  gas 
and  solid  density;  thus,  gas  phase  effects  were  ignored  in  this  study.  In  this  limit,  the  evolution  of 
mass,  momentum,  energy,  solid  volume  fraction,  and  no-load  (inelastic)  solid  volume  fraction  were 
described  by: 


9q  9g 

dt  dx  dy 


a) 
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i  T 
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otherwise. 


In  these  equations,  p  —  (j)ps  and  P  —  (j)Ps  are  the  effective  density  and  pressure  of  the  granular 
solid,  where  cj)  is  the  solid  volume  fraction  and  Ps  is  the  solid  pressure;  E  =  es  +  B  +  ( u 2  +  v2)/2  is 
the  mass-specific  total  energy  of  the  granular  solid,  where  es  is  the  solid  internal  energy  and  B  is 
the  compaction  potential  energy;  (f)  is  the  inelastic  component  of  the  solid  volume  fraction;  (3  is  a 
configurational  stress  that  resists  compaction;  and  /  is  a  yield  function  that  establishes  the  onset  of 
inelastic  volumetric  deformation.  The  parameters  p  and  p  establish  the  rates  of  stress  equilibration 
and  relaxation  to  the  yield  surface,  respectively. 

Of  particular  importance  in  this  analysis  is  the  heating  of  solid  explosive  because  of  its  relevance 
to  combustion.  For  low  pressure  waves,  the  dominant  heating  mechanism  is  volumetric  deformation 
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by  compaction,  where  the  mass-specific  compaction  induced  dissipative  heating  rate  is  given  by: 

wc  = 

{Pa-P)d4  /?# 
p  dt  p  dt 

Here,  d(*)/dt  =  d{9)/dt  +  v  •  V(«)  is  the  Lagrangian  derivative.  The  first  term  accounts  for  rate- 
dependent  dissipation  (it  vanishes  for  slow  compaction  because  Ps  «  /?),  whereas  the  second  term 
accounts  for  dissipation  due  to  inelastic  deformation.  It  is  noted  that  the  inclusion  of  <j6  in  the  theory 
enables  split  wave  structures  to  evolve  that  are  analogous  to  elastic-plastic  waves  in  dynamically 
loaded  homogeneous  solids  [4]. 

The  governing  equations,  expressed  in  terms  of  generalized  curvilinear  coordinates,  were  numer¬ 
ically  integrated  using  a  high-resolution,  shock-capturing  technique.  The  technique  is  nominally 
second-order  accurate  in  space,  and  third-order  accurate  in  time. 

2.2  Meso-Scale  Modeling  and  Computations 

Objectives.  Specific  objectives  of  the  meso-scale  modeling  and  computational  effort  included  the 
following: 

1.  Analysis  of  Micro- Cluster  Impact  with  a  Planar  Wall.  (See  Refs.  [7,  8,  9]  for  a  comprehensive 
discussion.)  The  objective  of  this  study  was  to  computationally  examine  inert,  plane  strain 
(2-D),  impact  of  a  small,  close,  randomly  packed,  and  numerically  well- resolved  micro-particle 
cluster  («  25  particles)  with  a  planar  boundary,  as  illustrated  in  Fig.  3(a).  Emphasis  was 
placed  on  characterizing  the  variation  in  temporal  and  spatial  partitioning  of  energy  within 
the  material  with  impact  angle  0,  and  on  identifying  the  relative  importance  of  plastic  and 
friction  work  as  potential  hot-spot  formation  mechanisms.  These  high  velocity  micro-clusters 
might  originate  in  practice,  for  example,  from  the  break-up  and  dispersal  of  energetic  solids 
by  strong  shocks.  This  analysis  is  important  for  the  following  reasons.  First,  the  impact 
angle  will  always  vary  for  loading  events  in  practice.  Second,  bulk  experiments  indicate 
that  combustion  of  energetic  solids  often  first  occurs  near  the  impactor  surface,  implying  that 
particle  heating  is  largest  at  that  location.  Third,  uncertainty  exists  about  the  coupled  effects 
of  particle-wall  and  particle-particle  frictional  interactions  on  stress  and  energy  states  within 
dynamically  loaded  materials  and  how  they  can  facilitate  reactive  hot-spot  formation.  Last, 
the  development  of  suitable  multi-dimensional  boundary  conditions  for  use  with  bulk  scale 
models  requires  a  firm  understanding  of  particle  scale  physics  occurring  in  the  immediate 
vicinity  of  solid  boundaries.  This  problem  provides  a  rational  foundation  for  additional  study 
on  the  influence  of  far- field  particle  interactions  (i.e.,  particle  self-confinement),  phase  change, 
and  combustion  on  the  3-D  near  wall  impact  response. 

2.  Analysis  of  Uniaxial  Deformation  Wave  Structure.  (See  Refs.  [8,  10]  for  a  comprehensive 
discussion.)  The  objective  of  this  study  was  computationally  examine  meso-scale  energy 
partitioning  within  uniaxial,  quasi-steady  deformation  waves  due  to  constant  speed  piston 
impact;  the  computational  domain  is  illustrated  in  Fig.  3(b).  Particular  emphasis  was  placed 
on  characterizing  the  evolution  of  hot-spot  mass  fraction  within  the  wave  structure  and  its 
dependence  on  impact  speed,  and  to  identify  the  relative  importance  of  plastic  and  friction 
work  as  hot-spot  mechanisms.  Quasi-ID  macro-scale  predictions  were  obtained  and  analyzed 
by  averaging  the  meso-scale  response  over  suitable  material  volumes.  A  further  objective  of 
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Figure  3:  Computational  configurations  used  for  the  meso-scale  analysis:  (a)  micro-cluster  impact; 
(b)  uniaxial  deformation  waves;  and  (c)  deformation  wave  interactions  with  non-planar  boundaries. 
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this  study  was  to  characterize  the  dependence  of  deformation  wave  structure  and  hot-spot  for¬ 
mation  on  constitutive  parameters  and  material  meso-structure.  Key  constitutive  parameters 
considered  included  the  plastic  viscosity,  which  establishes  the  time  scale  for  rate-dependent 
relaxation  to  the  material  yield  surface,  and  the  friction  coefficient,  which  establishes  the 
magnitude  of  the  surface  traction  at  inter-particle  contact  surfaces.  This  study  is  important 
because  both  of  these  parameters,  which  are  difficult  to  experimentally  determine  under  im¬ 
pact  conditions,  can  significantly  affect  the  particle  and  macro-scale  mechanics.  Also,  most 
energetic  materials  used  in  practice  have  a  non-uniform  particle  size  distribution  so  that  high 
volumetric  packing  densities  can  be  obtained.  It  may  also  be  possible  to  promote  a  desired 
effect  (e.g.,  impact  sensitivity  and  combustion  rate)  by  carefully  controlling  the  initial  particle 
size  distribution  of  the  material;  this  possibility  requires  a  reasonable  understanding  of  meso- 
scale  physics.  Therefore,  simulations  were  also  performed  for  materials  having  different  initial 
particle  size  distributions  to  demonstrate  their  effect  on  both  deformation  wave  structure  and 
hot-spot  formation. 

3.  Analysis  of  Deformation  Wave  Interactions  with  Non-Planar  Boundaries:  Comparison  with 
Macro-Scale  Predictions.  (See  Ref.  [5]  for  a  comprehensive  discussion.)  The  objective  of  this 
study  was  to  computationally  examine  the  2-D  inert  heating  response  of  granular  explosive 
due  to  interaction  of  an  initially  planar  compaction  wave  with  a  semi-circular  rigid  anvil;  the 
initial  computational  domain  is  illustrated  in  Fig.  3(c).  The  meso-scale  model  combined  con¬ 
servation  principles  with  an  elastic- viscoplastic  constitutive  theory  to  predict  contact  induced 
thermomechanical  fields  within  particles.  Meso-scale  field  predictions  were  locally  averaged 
and  compared  to  those  given  by  the  macro-scale  model.  Emphasis  was  placed  on  character¬ 
izing  the  evolution  of  spatial  wave  structure,  compaction  induced  heating  rate,  and  hot-spot 
mass  in  the  vicinity  of  the  anvil  surface,  and  their  dependence  on  initial  wave  strength.  This 
analysis  is  important  because  the  interaction  between  deformation  waves  in  granular  explo¬ 
sive  and  non-planar  macro-scale  boundaries,  which  commonly  occur  in  practice,  may  locally 
trigger  prompt  vigorous  combustion  due  to  high  dissipative  heating  rates,  as  mentioned.  A 
long  term  objective  of  this  work  is  to  correlate  meso-scale  hot-spot  predictions  with  macro¬ 
scale  quantities,  such  as  porosity,  pressure,  and  volumetric  deformation  rate,  for  use  with  a 
macro-scale  combustion  theory. 

Research  Method.  (A  comprehensive  discussion  of  the  numerical  modeling  technique,  and  its  ver¬ 
ification  and  validation  for  granular  HMX,  is  given  in  Refs.  [8,  9].)  Multi-body  contact  induced 
by  dynamic  compaction  was  posed  as  a  coupled  initial-boundary-value  problem  (IBVP)  for  the 
displacement  field  u  and  temperature  field  T  within  particles.  These  fields  were  described  by 
momentum  and  energy  evolution  equations  within  each  particle: 

pi i  =  V  •  er,  (3) 

pcvT  =  -V  •  q  +  pr.  (4) 

Here,  p  is  the  local  mass  density,  a  is  the  Cauchy  stress  tensor,  r  is  the  deformation  induced 
heating,  q  is  the  heat  flux,  cv  is  the  specific  heat  at  constant  volume,  and  V=  0(-)/dx  is  the  spatial 
gradient  operator.  Body  forces  were  ignored  as  they  are  inconsequential  compared  to  impact 
induced  deformation  forces.  All  particles  were  initially  stationary,  stress  free,  and  at  a  uniform 
ambient  temperature.  Only  contact  boundary  conditions  were  imposed  on  the  displacement  and 
temperature  fields  of  the  particles: 

<j  •  n  =  tc  onT  V  £, 

(— fcrVT)  •  n  =  qj  +  qc  onT  V  £, 
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(5) 

(6) 


where  T  is  the  particle  boundary,  n  is  the  unit  normal  to  T,  tc  is  the  contact  traction,  kx  is  the 
thermal  conductivity,  qj  is  the  heat  flux  due  to  frictional  heating,  qc  is  the  heat  flux  necessary  to 
impose  ideal  thermal  contact,  and  t  is  time.  The  system  of  equations  was  closed  by  prescribing  a 
constitutive  theory  for  the  stress  response.  To  this  end,  a  hyperelastic,  multiplicative,  finite  strain 
constitutive  theory  was  used  to  model  stress-strain  behavior.  A  Perzyna  type  over-stress  model, 
coupled  with  an  associative  flow  rule  and  a  Von-Mises  type  yield  criterion  with  isotropic  hardening, 
was  used  to  prescribe  the  evolution  of  inelastic  strain.  Material  properties  used  in  this  study  were 
chosen  to  be  representative  of  the  secondary  high  explosive  HMX. 

Simulations  were  performed  using  a  combined  finite  and  discrete  element  method  that  is  well- 
suited  for  problems  involving  heterogeneity.  This  combined  method  used  the  finite-element  method 
(FEM),  coupled  with  a  radial  return  stress  update  algorithm,  to  numerically  integrate  the  time- 
dependent,  2-D  conservation  principles  and  viscoplastic  flow  rule  governing  deformation  of  indi¬ 
vidual  particles,  and  used  the  discrete-element  method  (DEM)  to  account  for  interactions  between 
particles.  The  DEM  was  based  on  a  distributed,  conservative  potential  based  penalty  method 
whereby  the  normal  contact  traction  between  particles  was  estimated  by  penalizing  their  penetra¬ 
tion,  and  frictional  tractions  were  estimated  using  a  penalty  regularized  Amonton’s  Coulomb  law. 
Particles  were  discretized  using  constant  strain,  triangular  finite  elements,  where  each  particle  con¬ 
sisted  of  400-1200  elements  for  all  simulations  performed  in  this  work.  A  temporally  second-order 
accurate,  explicit  numerical  technique  was  used  to  integrate  the  finite  element  equations  for  nodal 
displacements  and  temperatures. 

3  Key  Results 

The  following  summary  highlights  key  results  for  the  objectives  identified  in  Section  2. 

3.1  Macro-Scale  Modeling  and  Computations 

1.  Analysis  of  Deformation  Wave  Interactions  with  Planar  Boundaries.  (See  Ref.  [6])  The  inert 
impact  response  of  initially  planar  deformation  waves  in  granular  HMX  with  rigid  planar 
boundaries  was  computationally  examined  [Fig.  2(a)].  Emphasis  was  placed  on  characterizing 
how  initial  wave  strength  and  wave  deflection  angle  6  affect  interaction  structure  and  explosive 
dissipative  heating.  To  this  end,  the  steady  1-D  equations  [obtained  from  Eq.  2],  valid  in 
a  wave-attached  frame,  were  numerically  integrated  as  an  Initial- Value-Problem  (IVP)  to 
obtain  the  spatial  wave  structure  corresponding  to  a  specified  impact  speed  Up\  the  resulting 
structure  was  then  interpolated  onto  the  2-D  computational  grid.  The  initial  radius  of  HMX 
particles  was  taken  to  be  r  —  100  /im,  and  the  initial  solid  volume  fraction  was  taken  to 
be  </>o  =  0.85.  Initial  wave  strengths  corresponding  to  Up  —  15,  120  and  300  m/s  were 
investigated  because  of  differences  in  their  spatial  structure:  the  structure  for  Up  =  15  m/s 
possesses  only  a  viscoelastic  region,  whereas  the  structures  for  Up  =  120  and  300  m/s  possess 
a  viscoelastic  precursor  followed  by  a  viscoplastic  region  (see  Ref.  [4]  for  a  comprehensive 
discussion  of  planar  deformation  wave  structures  in  granular  explosive).  The  wave  deflection 
angle  was  varied  between  0°  <  6  <  90°. 

Dispersed  wave-boundary  interaction  structures  were  predicted  that  are  analogous  to  those 
given  by  gas  dynamics  for  the  interaction  of  shocks  with  planar  boundaries.  As  illustrated 
in  Fig.  4  for  Up  =  120  m/s,  interaction  structures  generally  consisted  of  three  intersecting 
waves:  an  incident  wave,  a  stem  wave  near  the  boundary  surface  (analogous  to  a  Mach 
stem),  and  a  reflected  wave.  The  relative  strengths  of  these  waves  depend  on  Up  and  6.  For 
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Figure  4:  Predicted  mass-specific  dissipative  heating  rate  contours  for  the  interaction  of  an  initially 
planar  deformation  wave  with  a  rigid  planar  boundary  for  Up  =  120  m/s:  (a)  9  =  30°,  (b)  9  =  45°, 
(c)  9  —  60°,  and  (d)  9  —  75°.  The  spatial  contours  shown  do  not  necessarily  correspond  to  the  time 
of  maximum  heating. 


values  of  Up  considered  in  this  study,  these  waves  propagate  at  subsonic  speeds  relative  to 
the  ambient  pure  phase  explosive,  and  have  finite  thicknesses  due  to  deformation  induced 
dissipation  (unlike  gas  dynamic  waves).  The  wave  speeds  increase,  and  thicknesses  decrease, 
with  increasing  Up.  Borrowing  terminology  from  gas  dynamics,  Von  Neumann  reflections  are 
predicted  for  approximately  0°  <  9  <  30°,  Mach  reflections  are  predicted  for  approximately 
30°  <  9  <  60°,  and  regular  reflections  are  predicted  for  approximately  9  >  60°.  Deflection 
angles  which  define  these  different  flow  regimes  are  dependent  on  Up. 

The  dissipative  heating  rate  [given  by  Eq.  2]  is  particularly  important  because  of  its  relevance 
to  combustion  initiation;  it  is  predicted  to  increase  above  that  of  the  incident  wave  along 
the  surface  of  the  planar  boundary  due  to  the  interaction  process.  For  a  given  initial  wave 
strength  and  material  porosity,  the  peak  heating  rate  is  associated  with  Mach  reflections,  as 
illustrated  in  Fig.  4(c)  for  Up  =  120  m/s  and  9  =  60°.  Peak  heating  rates  of  wc  =  6.5,  13,  25, 
40,  16,  and  10  GW/kg  were  predicted  for  9  —  15°,  30°,  45°,  60°,  75°  and  90°,  respectively. 
The  peak  heating  rate  is  lower  for  9  =  75°  than  for  60°  because  the  stem  region  has  weakened 
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(a)  (b) 

Figure  5:  Predicted  spatial  variation  in  (a)  pressure  and  (b)  mass-specific  dissipative  heating  rate 
along  the  lower  boundary  of  the  computational  domain  for  Up  =  120  m/s  and  9  =  60°.  The  tip  of 
the  planar  boundary  is  located  at  x  —  10  mm.  Dash  curves  denote  reflected  waves. 


as  the  Mach  reflection  transitioned  to  a  regular  reflection  with  increasing  9.  Importantly, 
peak  heating  rates  were  predicted  to  occur  along  the  boundary  surface  near  locations  of  peak 
pressure  for  domains  investigated  in  this  study,  as  illustrated  in  Fig.  5  for  9  —  60°.  It  is  also 
noteworthy  that  reflected  waves  generated  by  wave-boundary  interactions,  which  propagate 
back  through  material  pre-compacted  by  the  incident  wave,  induce  little  additional  dissipative 
heating  which  is  indicative  of  shock  desensitization.  Predictions  for  other  piston  speeds  were 
qualitatively  similar  to  those  shown  here,  but  quantitatively  differed  in  terms  of  heating  rates 
and  wave  structure. 

Dissipative  heating  rates  due  to  rate-dependent  and  inelastic  deformation  [given  by  and 
e^  in  Eq.  (2),  respectively]  were  also  characterized,  which  depend  on  both  Up  and  9.  For 
Up  =  15  m/s,  dissipative  heating  by  the  incident  wave  was  solely  due  to  rate-dependent  (vis¬ 
coelastic)  deformation.  However,  Mach  reflections  were  shown  to  induce  additional  heating  by 
inelastic  deformation  along  the  boundary  surface  for  low  strength  waves  which  exceeded  that 
due  to  rate-dependent  deformation  for  large  9.  For  Up  =  120  and  300  m/s,  incident  waves  in¬ 
duced  significant  dissipative  heating  due  to  rate-dependent  deformation,  whereas  dissipative 
heating  due  to  both  rate-dependent  and  inelastic  deformation  were  shown  to  be  compara¬ 
ble  for  large  9.  These  results  indicate  that  wave-boundary  interactions  can  induce  complex 
macro-scale  dissipative  phenomena  that  may  separately,  or  in  concert,  trigger  combustion  of 
energetic  solids. 

2.  Analysis  of  Deformation  Wave  Interactions  with  Non-Planar  Boundaries.  (See  Refs.  [5,  6].) 
Simulations  were  performed  to  characterize  macro-scale  dissipative  heating  induced  by  in¬ 
teraction  of  initially  planar  deformation  waves  in  granular  HMX  with  a  rigid  circular  anvil 
boundary.  The  initial  radius  of  HMX  particles  was  r  —  70  /im,  and  the  initial  solid  volume 
fraction  was  </>o  =  0.835.  The  anvil  radius  was  taken  to  be  R  —  4  mm.  As  done  for  the  planar 
boundary  analysis,  the  initial  wave  structure  was  established  by  numerically  integrating  the 
appropriate  1-D  IVP  with  Up  specified,  and  interpolating  the  result  onto  the  2-D  domain. 

Figure  6  illustrates  the  predicted  spatial  variation  in  macro-scale  dissipative  heating  rate 
[defined  by  Eq.  (2)]  for  incident  waves  corresponding  to  Up  =  150,  300,  and  450  m/s.  Again, 
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Figure  6:  Predicted  mass-specific  dissipative  heating  rate  contours  for  the  interaction  of  an  initially 
planar  deformation  wave  with  a  rigid  circular  boundary  for  Up  =  150,  300,  and  450  m/s.  Both 
early  time  (left)  and  late  time  (right)  contours  are  shown. 
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interaction  structures  generally  consisted  of  an  incident  wave,  a  stem  wave  near  the  boundary 
surface,  and  a  reflected  wave.  Dissipative  heating  rates  increased  with  increasing  Up,  with 
peak  heating  occurring  within  the  stem  region  along  the  boundary  surface,  whereas  wave 
thicknesses  decreased.  Shown  in  Fig.  7  is  the  predicted  spatial  variation  in  heating  rate  at 
various  times  along  the  bottom  of  the  computational  domain  for  Up  —  300  m/s,  where  the 
stagnation  point  of  the  circular  anvil  is  located  at  x  =  8  mm;  results  were  qualitatively  similar 
for  all  incident  wave  strengths.  Two  observations  are  noteworthy.  First,  intense  heating  was 
induced  by  the  stem  wave  as  it  accelerates  along  the  front  surface  of  the  anvil  and  compacts 
virgin  material.  A  peak  heating  rate  of  wc  ~  800  GW/kg  was  predicted  to  occur  at  a  location 
along  the  anvil  surface  {x  ~  9.5  mm)  that  is  removed  from  the  high  pressure  stagnation 
region.  The  peak  heating  rate  was  nearly  fives  times  larger  than  that  of  the  incident  wave, 
and  would  likely  trigger  prompt  and  vigorous  combustion  of  reactive  material.  The  magnitude 
and  location  of  peak  heating  rate  was  dependent  on  anvil  size,  but  typically  occurred  near 
the  mid-point  of  the  anvil  surface.  Second,  the  reflected  wave  induced  low  heating  rates 
because  most  porosity  had  already  been  eliminated  by  the  incident  wave,  as  indicated  by  the 
dashed  history  curves  shown  in  Fig.  7.  In  this  respect,  the  pre-compacted  material  had  been 
desensitized  to  subsequent  wave  induced  deformation.  Though  not  shown  here,  peak  heating 
rates  of  wc  «  34  and  5000  GW/kg  were  predicted  for  Up  —  150  and  450  m/s,  respectively. 

It  should  be  re-emphasized  that  dissipative  heating  induced  by  wave  interactions  with  planar 
and  circular  boundaries  qualitatively  differed.  Peak  dissipative  heating  rates  and  solid  pres¬ 
sures  induced  by  wave  interactions  with  planar  boundaries  occurred  at  similar  locations  along 
the  boundary  surface,  whereas  peak  dissipative  heating  rates  induced  by  interactions  with  cir¬ 
cular  boundaries  occurred  at  locations  along  the  boundary  surface  that  were  well-removed 
from  the  high  pressure  stagnation  region.  These  predictions  are  consistent  with  the  exper¬ 
imental  observations  of  Ref.  [11]  which  showed  that  interaction  between  an  initially  planar 
shock  in  the  explosive  PBXW-128  and  a  hemispherical  glass  anvil  initiated  vigorous  com¬ 
bustion  near  the  explosive-anvil  boundary  away  from  the  high-pressure  stagnation  region, 
suggesting  that  maximum  heating  rate  occurred  at  that  location.  Moreover,  because  this 
analysis  ignored  the  effects  of  bulk  shear,  it  is  likely  that  combustion  initiation  may  be  pri¬ 
marily  due  to  volumetric  deformation,  with  bulk  shear  being  of  secondary  importance.  These 
predictions  collectively  indicate  that  dissipative  heating  rates  induced  by  wave-boundary  in¬ 
teractions  are  dependent  on  boundary  geometry,  and  that  pressure-dependent  models  used 
to  describe  combustion  of  energetic  solids  may  provide  misleading  results  for  such  problems. 

3.2  Meso-Scale  Modeling  and  Computations 

1.  Analysis  of  Micro- Cluster  Impact  with  a  Planar  Wall.  (See  Refs.  [7,  8,  9].)  Impact  of  an 
initially  stress  free,  close-packed,  2D  micro-particle  cluster  (25  particles  of  radii  50  fim)  and 
a  planar,  rigid  wall  was  computationally  examined  using  a  thermoelastic-viscoplastic  and 
stick-slip  friction  constitutive  theory.  The  influence  of  impact  angle  and  cluster  packing 
configuration  on  the  spatial  and  temporal  partitioning  of  dissipated  energy  was  characterized. 
As  illustrated  in  Fig.  8  for  a  uniform  impact  speed  300  m/s,  plastic  work  was  shown  to  rapidly 
increase  near  the  wall  immediately  following  impact  before  diffusing  outward  due  to  stress 
transmission  between  particles,  whereas  friction  work  near  the  wall  increased  sharply  after  a 
very  brief  initial  period  of  stick  («  75  ns  for  normal  impact).  The  initial  duration  of  stick  was 
largest  for  normal  impact  and  decreased  with  increase  in  impact  angle.  The  maximum  average 
plastic  work  was  approximately  4  kJ  per  unit  cluster  mass  for  an  impact  angle  of  <fr  =  0° 
(corresponding  to  normal  impact),  which  was  approximately  two  times  larger  than  that  for 


12 


Figure  7:  Predicted  spatial  variation  in  (a)  pressure  and  (b)  mass-specific  dissipative  heating  rate 
along  the  lower  boundary  of  the  computational  domain  for  Up  =  300  m/s.  The  stagnation  point  of 
the  circular  boundary  is  located  at  x  =  8  mm.  Dash  curves  denote  reflected  waves. 

(j)  —  60°.  The  maximum  average  friction  work  was  approximately  60  kJ  per  unit  cluster 
mass  for  <f>c  =  60°,  which  was  nine  times  larger  than  that  for  0  =  0°.  The  critical  impact 
angle  for  maximum  friction  work  was  shown  to  mildly  depend  on  the  initially  random  cluster 
configuration,  though  it  is  anticipated  that  little  variation  in  this  angle  would  result  for  larger 
particle  ensembles.  Friction  was  shown  to  significantly  affect  the  final  cluster  kinetic  energy, 
but  to  minimally  affect  elastic  potential  energy  and  plastic  work.  Also,  particle-wall  friction 
was  considerably  larger  than  inter-particle  friction  due  to  both  the  rigid  wall  assumption  and 
the  repeated  loading  of  particles  adjacent  to  the  wall  by  surrounding  particles.  Predictions 
from  the  energy  analysis  indicated  that  large  temperature  rises  (>  900  K)  occur  within  the 
cluster,  even  for  normal  impact,  with  the  maximum  temperature  (~  4400  K)  occurring  for 
4>  =  80°.  However,  the  mass  heated  to  such  elevated  temperatures  represented  less  than  .01% 
of  the  cluster  mass  with  the  most  mass  («  98%)  heated  to  temperatures  rises  <  200  K. 
Localization  of  dissipated  energy  indicated  that  average  temperature  rises  due  to  friction 
work  exceeded  those  due  to  plastic  work,  even  for  normal  impact,  and  far  exceeded  those 
due  to  plastic  work  for  highly  oblique  impact.  Average  frictionally  induced  temperature  rises 
near  900  K  were  predicted  for  (f>  —  0°  increasing  to  4400  K  for  </>  =  80°.  Impact  simulations 
for  a  deformable  wall  having  thermomechanical  properties  that  are  similar  to  the  cluster 
were  only  marginally  different  than  the  rigid  wall  results.  These  predictions  collectively 
highlight  the  importance  of  properly  describing  friction  as  a  local  heating  mechanism  that 
may  induce  combustion  of  energetic  particle  clusters.  This  study  has  provided  a  foundation 
for  systematically  characterizing  how  other  potentially  important  phenomena,  such  as  wall 
deformation,  particle  phase  change,  and  far-held  particle  interactions,  can  influence  the  near 
wall  impact  energetics  of  heterogeneous  solids. 

2.  Analysis  of  Uniaxial  Deformation  Wave  Structure.  (See  Refs.  [8,  10].)  Impact  induced  ther¬ 
momechanical  interactions  between  deformable  particles  contained  within  heterogeneous  en¬ 
ergetic  solids  were  computationally  examined  to  characterize  spatial  deformation  wave  struc¬ 
ture,  and  its  dependence  on  impact  speed  (vp  =  50-500  m/s),  and  to  identify  the  relative 
importance  of  friction  and  plastic  work  as  potential  hot-spot  formation  mechanisms.  Bulk 
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Figure  8:  Predictions  for  the  temporal  and  spatial  variation  in  average  specific  (a)  plastic  work  and 
(b)  friction  work.  Here,  x  is  normal  distance  from  the  wall. 


wave  structures  having  a  mostly  elastic  precursor  region,  followed  by  a  mostly  plastic  region 
in  which  the  porosity  is  significantly  reduced,  were  predicted  for  all  cases  considered  in  this 
study,  though  the  lengths  of  these  regions  varied  with  piston  speed,  as  illustrated  in  Fig.  9. 
For  low  piston  speeds  (vp  <  100  m/s),  the  wave  structure  was  predicted  to  be  thick  and  dis¬ 
persed  (^  30  average  particle  diameters),  where  compaction  was  largely  due  to  rigid  particle 
rearrangement  rather  than  plastic  void  collapse.  The  particle  ensemble  was  fully  compacted 
for  vp  >  400  m/s  due  to  significant  particle  plastic  deformation  associated  with  void  collapse. 
For  these  speeds,  the  wave  structure  consisted  of  a  very  thin  precursor  zone  (^  1  average 
particle  diameter)  followed  by  a  slightly  thicker  plastic  zone  (~  3  average  particle  diameters). 
For  piston  speeds  approximately  between  100  m/s  and  400  m/s,  local  spatial  variations  in 
porosity  within  the  meso-structure,  which  depend  on  the  initial  particle  size  distribution  and 
packing  configuration,  induced  large  amplitude  spatial  fluctuations  in  thermomechanical  fields 
behind  the  wave;  the  fluctuation  magnitudes  decreased  for  higher  impact  speeds  as  material 
strength  effects  become  less  consequential.  However,  it  is  noted  that  intra-particle  defects 
(e.g.,  voids  and  cracks),  which  were  not  modeled  in  this  study,  may  induce  high  frequency, 
large  amplitude  fluctuations  for  high  impact  speeds  as  the  bulk  wave  thickness  decreases. 

Plastic  work  was  shown  to  be  the  dominant  energy  dissipation  mechanism  in  that  it  constitutes 
a  substantial  fraction  of  total  system  energy  (25-35%  for  vp  =  50-500  m/s),  whereas  friction 
work  constitutes  a  smaller  fraction  (2-1%  for  vp  =  50-500  m/s).  The  maximum  plastic  and 
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Figure  9:  Predicted  finite  element  contours  of  Von  Mises  stress  re  within  particles  at  2.75  fi  s  for 
vp  =  50,  250  and  500  m/s.  Values  shown  in  the  color-bars  are  in  GPa. 


friction  work  occurred  close  to  the  piston  boundary  due  to  material  and  geometry  mismatches 
between  the  rigid,  planar  piston  and  the  initially  circular  particles.  Close  to  the  piston, 
friction  work  increased  from  0.2  J/m  to  1.3  J/m  for  vp  —  50  to  500  m/s,  whereas  plastic  work 
increased  from  0.3  J/m  to  45  J/m.  Predictions  for  mass-specific  friction  and  plastic  work 
highlighted  two  important  results.  First,  as  illustrated  in  Fig.  10,  mass-specific  friction  work 
dominated  plastic  work  for  vp  =  50  m/s,  whereas  plastic  work  dominated  friction  work  for 
vp  —  500  m/s.  This  result  suggests  that  friction  work  may  be  more  important  than  plastic 
work  in  producing  hot-spots  for  low  pressure  impact.  Second,  friction  work  is  sensitive  to 
details  of  the  meso-structure,  particularly  for  low  pressure  impact.  Predicted  maximum  hot¬ 
spot  temperatures  far  exceeded  their  average  values,  even  for  vp  =  500  m/s  where  a  maximum 
average  temperature  of  only  368  K  was  predicted;  this  average  value  which  is  well  below  the 
approximate  ignition  threshold  for  HMX  (Tig  ~  600  K).  Maximum  hot-spot  temperatures 
near  600  K  were  predicted  for  vp  =  50  m/s,  which  increased  to  near  1100  K  and  1400  K  for 
vp  —  250  m/s  and  500  m/s.  Meso-scale  temperature  fields  indicated  that  hot-spots  induced 
by  both  plastic  and  friction  work  are  of  sub-particle  size.  Frictionally  induced  hot-spots 
occurred  along  contact  interfaces;  plastically  induced  hot-spots  occurred  within  particles  in 
the  vicinity  of  contact  surfaces,  and  they  were  substantially  larger  in  size  and  greater  in 
number  than  frictional  hot-spots. 
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Heated  mass  was  predicted  to  steadily  increase  with  distance  behind  the  wave  front  until  a 
quasi-steady  value  was  reached  at  the  end  of  the  compaction  zone,  with  minimal  variation  with 
position  occurring  beyond  that  location.  Minimal  heating  was  predicted  within  the  precursor 
region  of  the  compaction  zone  which  resulted  in  temperature  rises  ofA  T  <  5  K  for  all  cases 
considered.  As  illustrated  in  Fig.  11,  quasi-steady,  local  hot-spot  mass  fraction  distribution 
curves  behind  the  wave  consisted  of  three  parts  based  on  temperature  rise  intervalsA  T ^  and 
A Tf.  The  first  part  included  mass- fractions  ranging  from  10-1  and  10-3  that  were  heated 
through  a  temperature  rise  ofA  T5  «  30,  100,  and  200  K  for  vp  =  50,  250,  and  500  m/s. 
These  hot-spots  primarily  resulted  from  plastic  work  occurring  within  the  interior  of  particles. 
The  second  part  of  the  distribution  curves  included  mass-fractions  ranging  between  10-3  to 
10-4  that  were  heated  through  a  temperature  rise  betweenA  T &  andA  Ty ,  where  A  Tf  ~  70, 
200,  and  400  K  for  vp  —  50,  250,  and  500  m/s,  respectively.  These  hot-spots  resulted  from 
a  combination  of  plastic  and  friction  work  occurring  in  the  vicinity  of  contact  surfaces.  The 
third  part,  or  tail,  of  the  distribution  curves  consisted  of  small  mass-fractions  (m  =  10-5) 
that  experiencedA  T  >  A  Tf.  These  hot-spots,  which  resulted  from  friction  work,  varied 
sporadically  with  axial  position  behind  the  wave  due  to  spatially  nonuniform  heterogeneities 
within  the  meso-structure;  they  increased  in  number  with  piston  speed.  By  comparing  hot¬ 
spot  contours  obtained  using  the  predicted  temperature  field  to  those  estimated  using  an 
adiabatic  plasticity  analysis,  where  frictional  effects  and  thermal  conduction  were  ignored, 
revealed  that  only  highly  isolated  mass  (mass-fraction  10  5)  experienced  AT  >  25  K  due 
to  plasticity  for  vp  =  50  m/s.  For  this  speed,  almost  all  mass  heated  above  25  K  was  entirely 
due  to  frictional  dissipation.  For  vp  =  250  m/s  and  500  m/s,  almost  all  mass  heated  to 
AT  >  300  K  was  due  to  a  combination  of  plasticity  and  friction,  or  friction  alone.  These 
predictions  collectively  highlight  that  friction  work  is  the  more  important  hot-spot  mechanism 
for  low  speed  impact,  and  though  often  ignored,  is  also  significant  for  high  speed  impact. 

Sensitivity  studies  were  performed  to  examine  the  influence  of  friction  coefficient,  /i,  on  com¬ 
paction  wave  structure  by  varying  its  value  from  zero  to  the  baseline  value  of  0.25.  Predictions 
showed  that  the  average  porosity  behind  waves  marginally  decreased  with  decreasing  friction 
coefficient  due  to  enhanced  particle  mobility.  Increasing  the  friction  coefficient  promoted 
stick  behavior  between  particles  resulting  in  a  more  frictionally  rigid  ensemble.  Frictional 
stick  decreased  particle  mobility  at  higher  friction  coefficients  which  promoted  more  efficient 


Figure  10:  Predictions  for  the  partitioning  of  average  mass-specific  system  energy  at  2.75  /x  s  for 
(a)  vp  —  50  m/s  and  (b)  vp  =  500  m/s.  In  these  plots,  average  mass-specific  plastic  work  is  given 
by  wp ,  and  average  mass- specific  friction  work  is  given  by  Wf. 
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Figure  11:  Predicted  spatial  variation  of  hot-spot  mass- fraction  through  the  compaction  wave  at 
2.75  /x  s  for  vp  =  50  and  500  m/s.  Plots  to  the  left  give  the  average  rates  of  plastic  and  friction 
work.  The  plots  to  the  right  give  hot-spot  mass  fraction  at  locations  indicated  in  the  corresponding 
plots  to  the  left. 

stress  transmission  between  particles  and  faster  bulk  wave  speeds.  Predicted  pressure  and 
stress  fields  showed  minimal  sensitivity  to  changes  in  friction  coefficient.  Though  average  tem¬ 
peratures  were  predicted  to  marginally  increase  with  friction  coefficient,  maximum  hot-spot 
temperatures  were  shown  to  increase  substantially.  Peak  temperatures  of  360,  460  and  660  K 
were  predicted  for  frictionless  cases  corresponding  to  vp  =  50,  250  and  500  m/s,  respectively, 
whereas  temperatures  close  to  600,  1100  and  1400  K  were  predicted  for  cases  have  a  baseline 
value  of  fi  —  0.25.  Predictions  also  showed  that  frictionally  induced  hot-spot  mass-fractions, 
which  constitute  the  high  temperature  end  of  distributions,  increased  significantly  with  an 
increase  in  friction  coefficient,  as  illustrated  in  Fig.  12  for  vp  =  250  m/s. 

Sensitivity  studies  were  also  performed  to  examine  the  influence  of  plastic  viscosity,  7,  on 
compaction  wave  structure  by  varying  its  value  from  10  Pa-s  to  its  baseline  value  of  100  Pa-s. 
Lowering  viscosity  increased  inelastic  deformation  rates  and  resulted  in  lower  porosity  behind 
waves  and  marginally  lower  wave  speeds.  Predicted  pressure  fields  were  shown  to  be  insensi¬ 
tive  to  changes  in  viscosity.  Lower  viscosity  also  resulted  in  lower  average  Von  Mises  stresses 
re  but  higher  plastic  strains  ep.  Thus,  plastic  work,  given  by  the  product  of  te  and  ep,  was 
only  marginally  affected  by  the  viscosity.  Predicted  average  temperatures  were  insensitive 
to  viscosity,  though  peak  temperatures  increased  with  viscosity  for  vp  >  100  m/s.  Predic¬ 
tions  showed  that  viscosity  can  substantially  enhance  the  formation  of  frictionally  induced 
hot-spots,  as  illustrated  in  Fig.  13. 
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Figure  12:  Predicted  hot-spot  mass-fraction  contours  for  vp  =  250  m/s  at  2.75  /jl  s  for  (a)  ji  —  0.0 
and  (b)  \i  —  0.25.  Values  in  the  color-bar  represent  the  logarithm  of  mass-fraction.  The  superposed 
black  curves  represent  the  predicted  spatial  profile  for  the  average  rate  of  plastic  work. 


Finally,  sensitivity  studies  were  performed  to  examine  the  influence  of  initial  particle  size  dis¬ 
tribution  (PSD)  on  deformation  wave  structure.  To  this  end,  an  alternate  particle  ensemble 
was  generated  that  had  the  same  initial  porosity  as  the  baseline  ensemble,  but  a  wider  PSD. 
Predictions  for  the  mechanical  response  of  both  ensembles  showed  similar  qualitative  and 
quantitative  trends.  The  wider  PSD  resulted  in  marginally  lower  porosity  behind  waves,  and 
wave  speeds,  due  to  increased  mobility  of  smaller  particles.  For  a  given  impact  speed,  quasi¬ 
steady  wave  pressures  for  both  ensembles  were  similar,  but  the  stresses  within  smaller  particles 
exhibited  larger  spatial  fluctuations  than  those  within  larger  particles.  The  inelastic  defor¬ 
mation  and  average  temperature  responses  of  both  ensembles  slightly  differed,  with  smaller 
particles  exhibiting  higher  average  temperatures  than  larger  particles.  The  frictional  response 
was  shown  to  be  highly  sensitive  to  PSD  in  that  a  wider  PSD  enhanced  frictional  dissipation 
resulting  in  higher  heating  rates  and  temperatures  for  all  piston  speeds,  as  illustrated  in  Fig. 
14.  Though  smaller  particles  experienced  higher  average  temperatures,  predicted  hot-spot 
mass  fraction  distributions  for  each  particle  size  class  indicated  that  larger  particles  experi¬ 
enced  localized  heating  which  more  significantly  contributed  to  the  high  temperature  end  of 
the  distributions.  These  predictions  suggests  that  larger  particles  may  play  a  prominent  role 
in  combustion  initiation  of  energetic  solids  whereas  smaller  particles  may  be  important  for 
reaction  propagation. 

3.  Analysis  of  Deformation  Wave  Interactions  with  Non-Planar  Boundaries:  Comparison  with 
Macro-Scale  Predictions.  (See  Ref.  [5].)  Simulations  were  performed  to  computationally 
examine  multi-scale  interaction  of  initially  planar  deformation  waves  in  granular  HMX  with  a 
semi-circular  rigid  anvil.  The  analysis  focused  on  low  pressure  waves  resulting  from  100-450 
m/s  impact.  The  initial  solid  volume  fraction  of  the  material  was  taken  to  be  </>o  =  0.835,  and 
the  average  grain  size  was  ro  =  70  /am.  The  anvil  radius  was  taken  to  be  R  =  4  mm.  The 
domain  size  was  largely  chosen  based  on  computational  time  constraints,  and  to  maintain 
consistency  with  the  macro-scale  analysis;  a  typical  simulation  required  three  days  on  a 
parallel  computer  architecture. 

Figure  15  gives  the  average  mass-specific  dissipative  heating  rate  predicted  by  the  meso-scale 
model  for  Up  =  150,  300,  and  450  m/s.  Meso-scale  averages  were  locally  computed  over 
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Figure  13:  Predicted  hot-spot  mass-fraction  contours  for  vp  —  500  m/s  at  2.75  fi  s  for  (a)  7  =  1  x  101 
Pa-s  and  (b)  7  =  1  x  102  Pa-s.  Values  in  the  color-bar  represent  the  logarithm  of  mass-fraction. 
The  superposed  black  curves  represent  the  predicted  spatial  profile  for  the  average  rate  of  plastic 
work. 


a  square  averaging  area  of  500  x  500  fiw?  to  facilitate  comparison  to  macro-scale  model 
predictions.  Here,  the  dissipative  heating  rate  was  due  to  viscoplastic  deformation  of  grains 
during  pore-collapse.  For  all  values  of  Up,  average  meso-scale  predictions  qualitatively  agree 
with  macro-scale  predictions  in  that  peak  heating  rates  are  induced  by  the  stem  wave  at  a 
location  removed  from  the  stagnation  region.  Also,  pre-compacted  material  was  desensitized 
by  the  incident  wave.  However,  for  approximately  Up  >  300  m/s,  the  magnitude  of  heating 
rates  predicted  by  the  meso-scale  model  were  significantly  lower  than  those  of  the  macro¬ 
scale  model,  and  wave  thicknesses  were  larger.  For  these  cases,  the  averaging  area  width 
exceeded  the  wave  width  causing  the  average  dissipative  heating  rate  within  the  wave  to 
be  unresolved  at  the  macro-scale.  For  approximately  Up  <  200  m/s,  meso-  and  macro-scale 
predictions  agreed  well  both  qualitatively  and  quantitatively  based  on  the  averaging  area  used 
in  this  analysis.  Establishing  the  proper  choice  of  averaging  area  size,  and  its  dependence  on 
material  meso-structure  and  loading  conditions,  is  an  important  topic  of  ongoing  research. 

Meso-scale  simulations  also  provided  detailed  information  about  local  fluctuations  in  thermo¬ 
mechanical  fields  that  are  important  for  combustion  initiation.  Of  particular  importance  was 
the  evolution  of  hot-spot  mass  fraction  distribution.  Figure  16  gives  predicted  distributions 
for  averaging  areas  located  within  the  stagnation  and  maximum  heating  regions.  Here,  m 
is  the  mass  fraction  and  A  T  =  T  —  To,  where  To  is  the  ambient  temperature.  Though  the 
low  temperature  end  of  the  distributions  (AT  <  200  K)  were  similar  in  the  stagnation  and 
maximum  heating  rate  regions,  the  high  temperature  ends  departed  indicating  a  substan¬ 
tial  increase  in  hot-spot  mass  within  the  latter  region.  This  difference  in  hot-spot  mass  is 
significant,  and  will  likely  affect  the  explosive  combustion  response.  These  predictions  are 
again  consistent  with  the  experimental  observations  of  Ref.  [11].  The  coupled  dissipative 
heating-combustion  response  of  the  explosive  material  is  an  ongoing  topic  of  our  research. 
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Figure  14:  Predicted  hot-spot  mass-fraction  contours  for  vp  =  50  and  500  m/s  at  2.75  /a  s  for  config¬ 
uration  S-l  and  S-2.  Values  in  the  color-bar  represent  the  logarithm  of  mass- fraction.  Configuration 
S-l  represents  the  baseline  ensemble,  and  configuration  S-2  represents  the  alternate  ensemble  having 
a  wider  PSD. 
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Figure  15:  Meso-scale  predictions  for  the  average  inelastic  heating  rate  induced  by  the  interaction 
of  an  initially  planar  deformation  wave  and  a  semi-circular  macro-scale  boundary. 
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Figure  16:  Predicted  hot  spot  mass  fraction  distribution  in  the  stagnation  and  maximum  heating 
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